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When the process of a system in contact with a heat bath is described by classical Langevin 
equation, the method of stochastic energetics [K. Sekimoto, J. Phys. Soc. Jpn. 66 (1997) 1234] 
enables to derive the form of Helmholtz free energy and the dissipation function of the system. 
We prove that the irreversible heat Qi rr and the time lapse At of an isothermal process obey the 
complementarity relation, Q[ TI At > k^T 5 m in, where <S m i n depends on the initial and the final 
values of the control parameters, but it does not depend on the pathway between these values. 
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The relations of thermodynamics tell that, when we 
change quasi-statically the control parameters of a sys- 
tem in contact with a heat bath, the work W needed for 
the change is equal to the change of Helmholtz free en- 
ergy, AF, which consists of the reversible heat released 
to the heat bath, Q Tev , and the change of internal en- 
ergy, AE, as Q re v + AE = AF. When the change of the 
control parameter is not quasi-static, the needed work is 
more than the reversible one, i.e., W — AF > 0, and the 
released heat Q, which obeys the energy conservation 
law Q + AE = W, is larger than Q rev by the amount 
so-called irreversible heat Q m = Q — Q Tey = W — AF. 

In order to assess the released heat Q under a given 
protocol of control parameters, we need both a dynam- 
ical model of a system and a proper kinematical inter- 
pretation of the heat release from the system. One of 
the present authors© have introduced a method to ob- 
tain Q in the systems whose dynamics is described by 
Langevin equations. We will show below that, if applied 
to slow change of control parameters, this method, which 
we shall provisionally call stochastic energetics, enables 
to formulate reversible and irreversible thermodynamics 
of processes of the system in contact with a heat bath. 

First we recapitulate the main idea of stochastic en- 
ergetics in the case of a single heat bath. Let x — 
{xi, . . . ,x n } represent the state of the fluctuating sys- 
tem and let a — {a\, . . . , a r } be the parameters which 
control the system through the potential U(x;a). The 
Langevin equation is assumed as follows, 

dx . . dU , 
= -r- — +£(*)-— (x;a), (1) 



dt 



where F is a positive definite and symmetric friction ma- 
trix, which we assume to be constant, and is white 
and Gaussian random forces characterized by (£(t)) = 
and (£(*]*£(*')) = 2Tk B Td(t - t')B As far as we can 
regard (|l|) as a mechanical balance equation, the bal- 



ance equation for energy is obtained by making the scalar 
products of each term with dx{t) = x(t+ ^j)— x(t— 
along the realized trajectory, 



dU . , . , ,. 
— (x(t);a(t)) 



*(*)) ' dx(t) 

da(t)-dU(x{t);a{t)), (2) 



§R . dx+ 
dx ajj ^ da 



where we have used the identity dU 
da, which is valid for stochastic variables as far as the 
multiplication is defined in the Stratonovich sense.™ 

The left hand side is the heat released to the heat 
bath, dQ, since — (— T%(t) + £(t)) is the reaction force 
which the system exerts onto the heat bath. The energy 
conservation law implies that dQ + dU should be equal to 
the work dW done by the external agent to the system, 
hence 

dW = ^(x{t);a{t))-da(t). (3) 

Suppose that the control parameters a is changed (by 
some external agent) from a t at t = to <2f at t — At. 
The total work W done to the system in the course of a 
particular process x(t) (0 < t < At) is then, 



w = 



, dU , , , . . , 
dt—(x(t);a(t)) 



da(t) 
dt 



(4) 



That the work W is defined before we take ensemble 
average is both conceptual and practical advantage of 
stochastic energetics as compared with the master equa- 
tion approaches .Q'cP The ensemble average of W over 
possible realization of {£(<)}o<t<At is expressed as 



W={W) 



A/ 



dt 



J dxP{x,t)^-{x;a(t)) 



da(t) 
dt 



, (5) 



where P is the probability distribution function of x 
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which obeys the following Fokker-Planck equation^ 
dP 



dt 



(x,t) = -C FF (a(t))P(x,t) 

d tfdTI d 

-•I- (-(x M t ))+ k B T-)P(x,t).(, i 



Below we describe the main results, whose derivation 
will be given in the last part of the text. Let d(s) be 
a given protocol of the control parameter a with a unit 
time lapse (At = 1) satisfying d(0) = ctj and d(l) = Of. 
We assume, for convenience, the unit of time to be the 
characteristic time of equilibration of P(x, t) when a is 
fixed at a typical value of the protocol. Then for the 
stretched or slowed-down protocol a(t) = a(-^) with a 
large time lapse At, the following asymptotic formula for 
W holds, 



W = AF 



+ - 



k B T 



At 



ds 



t/ da(s) 
ds 



A(d( S )) 



da(s) 
ds 



+0((Aty 2 ), 



(7) 



where AF is the difference of the Helmholtz free energy, 
AF = F(a f )- F(di) with 



F(a) = —k B T log I j dx exp 



U(x; a) 
k B T 



const. 



(8) 

and A(a) is a positive definite n x n matrix defined by 



A(o) 



dx / dx' 



dP 



cq 



(x;a)g(x,x';a) 



tfdP 



cq 



da 



(x';a) 



(9) 



where P cq (x;a) = e -u{x-,a)/k B T j j Me -u(x-,a)/k B T ig 
the equilibrium distribution under a given parameter 
value, a, and the kernel g(x,x';a) is defined as the so- 
lution of the following equation 



- £Fp(a)[P eq (x;a)g(x,x'; a)} = 5(x - x') 



(10) 



The result (7) tells firstly that, in the limit of slow 
and smooth change of external parameters (At — > oo), 
the stochastic energetics reproduces correctly the ther- 
modynamic relation of quasi-static isothermal processes, 
W = AF. Secondly, for large but finite time lapse At, 
the irreversible heat Q m — W — AF behaves asymp- 
totically as ~ 4t, with the proportionality constant de- 
pending on the scaled protocol a(-^). The integral on 
the right hand side of (7) is analogous to an classical ac- 
tion of free particle with a 'mass' 2A being dependent 
on the 'position' a. The minimum value of this integral, 
which we denote by <S m i n (ai, a.f), should be realized by a 
certain 'classical' path, or the optimal (scaled) protocol 
d*(s) (0 < s < 1). We then come to the following com- 
plementarity relation which is correct asymptotically for 
A< -> oo, 



Q irr At > k B TS min (ai, a f ) 



(11) 



The above relation implies a sort of uncertainty relation 



that the estimation of the Helmholtz free energy function 
by the measurement of mechanical work W inevitably in- 
cludes a deviation Q m whose lower bound is controlled 
by the inverse of the measuring time At, at least in mod- 
erately slow measurements. 

Remark 1. The dissipation function $ of linear irre- 
versible thermodynamics is usually defined as A 3'" = 
2 $(»). In our method $ is obtained if we neglect the 
0((At) 2 ) terms in (7) and rewrite the integral there 
using a(t) instead of a(s); 



da(t) k B tfda(t)\ 



dt 



-* 2 1 dt ) 



A(a(t)) 



da(t) 
dt 



(12) 



Remark 2. In terms of the spectral representation of 
the Fokker-Planck operator, 



C FP (a) = 22 \m:a)X m (a)(m;a\ 



(13) 



with = Ao(a) < A x (a) < X 2 (a) < . . . and (m; a\n; a) 
5 mn , the kernel A(a) is formally expressed as 

1 



A(o) 



m— 1 



(k B T) 
(0;a 



t/dU 



V da 



\m;a) \ (m;a\^\0:a 
Am (a) da 



«rl0;a>-(i4) 



Here a appears merely as a parameter. 

Remark 3. The function g(x, x 1 ; a) is the Green's func- 
tion of a Hermitian operator; 

kBT iL ■ r ~ ipc ^ x; °) • a( - x > x '> °) = 5{ - x - x '^- 

(15) 

Especially when a; is a single variable, x, and T, a con- 
stant 7, the explicit form of g(x, x' ; a) is 

g(x, x'; a) = — Y — sgn(x - x') [ ^ (16) 
2fc B T J x , P eq (z;a) 

Remark 4- If the control parameters a are constrained 
to change along a single trajectory, say a(9), where 
< 8 < 1 is the contour parameter with d(0) = ai 
and d(l) = af, then Qj rr depends yet on how 9 depends 
on the scaled time s, as well as on the time lapse At. The 
integral in (7) then becomes J A(9(s))\^(s)\ ds, where 

A(0) = t (^) • A(o(0)) • ^p-. If we note that this 

integral becomes the simple action \j^(s)\ 2 ds by the 

transformation du(s) = A(9(s)) 1 ^ 2 d9(s), the minimum 

1 ^ 

of the integral is explicitly given as K(9) X ^ 2 d9 . Al- 
though not apparent, this result is invariant under the 
relabeling of the contour parameter. 

As an illustration, let us apply our formulae to the case 
of U(x; a) = 7}X 2 with the single variable x and the single 
control parameter a, with T a constant 7. In O(At ) 
order we recover the correct (configurational) free energy 

In this case F is purely of entropic origin since, from 
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cquipartition theorem, the internal energy is independent 
of the 'width' of the potential, ~ \ j\fa. 

In 0((Aty L ) order, we obtain the irreversible heat 



Q irr 



4Af Jo ~"a(s) 



ds- 



da(s) 



ds 



o((Aty 



(18) 



The integral can be minimized by the following scaled 
protocol; 

-2 



a(s)\ 



ptimum 



S 

at 



1 - s 



(19) 



and the asymptotic complementarity relation reads as 
follows, 

2 



r At > 



1 



1 



(20) 



Proof of (7) If we introduce the scaled time s = t/ At 
and the probability distribution with this argument, 
P(x, s; At) = P(x,sAt), the equations (5) and(6) be- 
comes, respectively, 



W 



da^Jfl . f dx^-(x; a(s))P(x,s; At) (21) 



ds 



da 



1 dP 

— — (a;, s; At) = -£ FP (a (s P (x, s; At). (22) 
At os 



P (1 \x,s) = P cq (x;a(s)) 
dP c 

dx'g(x, x'; a{s))—^{x'; a(s)) + X 



> (25) 



where \ lS to be determined from the normalization con- 
dition, J dxP^\x, s) — 0. After performing this, and 

noting the relation ^f-(x' ;a(s)) = (x' \ a(s)j\ ■ 

, the kernel A(d(s)) in (7) becomes as follows, 



A.(a) = -j^; J dx J dx' t (^-(x;a)^J P cq {x;a 



BP 
oa 



(26) 



where we have introduced an operator Vx (o) such that 
V^(a)^j(x) = J dx[S(x -x)- P cq {x; a)]rb{x) (27) 



for an arbitrary function of ip(x). Next the factors 
(Hj(a;;a)) P cq (x;a) in (26) can be rewritten as 

(jj^( x ': a )^j Pcq(x; a) = -k B T^^(x; a) 



P eq (x;a) J dx (^-(x:a)^j P cq (x;a). 



(28) 



Here the second term on the right hand side 
(oc P cq (x; a)) does not contribute to the integral of (26) 
since the following identity holds for an arbitrary func- 
tion, ip(x), 



J dxP eq (x;a)[V x -(a)iP(x)] = 0. 



(29) 



Now the integral in A(a) is of the form of 
J dx^ a ± {x;a)ip{x), and since the following identity 



f BP 
J dx-^(x;a)[P^(a)^(x)] 



BP C 

dx g Cq ( x i a )i } { x ) 



(30) 



For large At, we may solve (|22|) perturbatively in the 
form of 

P(x, s; At) = (x, s) + ^P (1) (x, s) + • • • . (23) 

In the lowest order p(°) obeys Cfp(o,(s)JP^ ' (x, s) = 0, 
and the normalization condition j dxP^ (x, s) = 1. 
Then, p(°\x,s) is the equilibrium distribution for a 
given parameter d(s), i.e., 

PM(x,s) = P eq (x;a(s)), (24) 

From (5) and the identity J dx§^(x; a)P cq (x; a) = 
|f (a), the O(At ) term of W becomes F(a f )- F( ai ). 

In the next order, (p2) becomes, — £fp(Q'(s))P^ 1 - ) (sc, s) = 
J|P eq (a;, d(s)). Using (!(]), the solution is given as 



holds for an arbitrary function ip(x), we come to the 
expression (9). 

In this letter we have considered the irreversible pro- 
cess near equilibrium states. Although it may not have 
been explicitly stated, the complementarity relation itself 
is derivable from phenomenological dissipation function. 
The method of stochastic energetics not only gives a con- 
crete form of the bound of the uncertainty between Q ilr 
and At, but it is also applicable to not slow processes 
such as thermal ratchetsEP or far from equilibrium non- 
steady statesH 1 How the irreversible heat does or does 
not depend on the models under the non-slow processes 
is a challenging problem in future.™ 
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